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ABSTRACT 

We present multi-point in situ observations of a complex sequence of coronal mass ejections (CMEs) 
which may serve as a benchmark event for numerical and empirical space weather prediction models. 
On 2010 August 1, instruments on various space missions (Solar Dynamics Observatory/ Solar and He- 
liospheric Observatory/Solar-TErrestrial-RElations-Observatory) monitored several CMEs originating 
within tens of degrees from solar disk center. We compare their imprints on four widely separated 
locations, spanning 120° in hcliospheric longitude, with radial distances from the Sun ranging from 
MESSENGER (0.38 AU) to Venus Express (VEX, at 0.72 AU) to Wind, ACE and ARTEMIS near 
Earth and STEREO-B close to 1 AU. Calculating shock and flux rope parameters at each location 
points to a non-spherical shape of the shock, and shows the global configuration of the interplanetary 
coronal mass ejections (ICMEs), which have interacted, but do not seem to have merged. VEX and 
STEREO-B observed similar magnetic flux ropes (MFRs), in contrast to structures at Wind. The 
geomagnetic storm was intense, reaching two minima in the Dst index (« —100 nT), caused by the 
sheath region behind the shock and one of two observed MFRs. MESSENGER received a glancing 
blow of the ICMEs, and the events missed STEREO-A entirely. The observations demonstrate how 
sympathetic solar eruptions may immerse at least 1/3 of the heliosphere in the ecliptic with their 
distinct plasma and magnetic field signatures. We also emphasize the difficulties in linking the local 
views derived from single-spacecraft observations to a consistent global picture, pointing to possible 
alterations from the classical picture of ICMEs. 

Subject headings: shock waves - solar-terrestrial relations - Sun: coronal mass ejections (CMEs); Sun: 
heliosphere 



1. INTRODUCTION 

Interplanetary coronal mass ejections (ICMEs) are 
massive expulsions of plasma and magnetic flux from 
the solar corona into interplanetary space and are 
thought to form enormous bent tubes, with radial scale 
sizes of tenths of an AU, when they are observed 
in situ by spac ecraft in the solar wind. With few 
exceptions (e.g.. iBurlaea et al.l ri981b iCane et al.l 119971: 
Bothmer fe Schwennlll998t ILiu et aljl2008al iMostl et al l 
2009alld: iRouillard et aiT 120101: iFarrugia et al.l 120111: 
Kilpua et al.l 120111 : iRuffenach et al.l 120121 ). these mea- 
surements have been largely restricted to single-point 
observations in spac e. Since their discov ery in the be- 
ginning of the 1980s (jBurlaga et al.l lT981). this situation 



has left researchers wondering about their global config- 
uration, because the extreme undersampling results in 
too many free parameters to unambiguously determine 
their shape, geometry and magnetic field topology on a 
large scale. The same is true for shock waves in the solar 
wind that are often driven by the magnetic str ucture in- 
side the ICME fe.g. JRichardson fe Cane1ll993D . where it 
is in general also unclear how to extrapolate, if possible 
at all, the global shock shape and structure from pa- 
rameters obtained only at a single point in space. This 
is important because the structure of th e shock deter- 
mines how par ticles are accelerated (e.g.. iReame si IT999I : 
ILiu et al.ll20lTI ). 

Furthermore, for testing real-time space weather pre- 
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diction models (e.g.. iSiscoe fc Schwenn|[2006h . it would 
be highly desirable to pinpoint the exact arrival times 
of an ICME at several locations in the inner helio- 
sphere. Using these known arrival times as bound- 
ary conditions, researchers can fine -tune numerical sim- 
ulations (e.g., Od strcil et al.l 120041) o r empirical meth- 
ods (e.g.. iGopalswamv et all l200ll Kim et al.l l2Q07t 
IVrsnak et al.l 12012ft . in order to improve the predic- 
tion of the ICMEs shock and magnetic ejecta arrival 
time at Earth. On a deeper level, we ultimately want 
to understand the physics of the ICME's motion in 
the solar wind, which was found to be dependent on 
their heliospheric e nvironment (e.g.. lTemmer et al.ll201 it 
iKilpua e~ al. 2012). Single-spacecraft in situ observa- 
tions can constra in ICME propagation models (e.g., 
iRollett et al.l [20 12). but it is often unclear if the prin- 
cipal direction of the ICME at 1 AU is east or west 
of the observer and how the shock surface and the flux 
rope "tube" look on a global scale, which is needed for 
forecasting the arrival time wit h geometrical models ap- 
proximating these s hapes (e.g.. [Howard fe Tappinl 120091: 
Lugaz et all 120091; [Mostl et al.l 120111 : iMostl fc Daviesl 
2012t iDavies et al 1120121 ). With the advent of the Solar- 



TErrestrial-R Elations-Observato ry (STEREO, launched 
in late 2006, Kaiser et al.l 120081) and the Solar Dynam- 
ics Observatory (SDO, launched in early 2010), a new 
era has begun where in situ multi-point observations of 
ICMEs have become available in connection with high- 
resolution imaging of the solar corona and multi-point 
imaging of the entire Sun-Earth line. 

Historically, because of a lack of systematic measure- 
ments, there have been very few in situ observations of 
ICMEs at large longitudinal or in-eclipti c separations, 
and those found were serendipito us (e.g., Burlag a et al.l 
[l98lt iBothrner fc Schwennl 11998ft . The STEREO mis- 
sion, consisting of two spacecraft increasingly leading and 
lagging the Earth at a rate of 22° per year, has made 
some observations of magnetic flux ropes (MFRs; inter- 
vals with smoothly rotating magnetic field) inside ICMEs 
at up to three points in space (including Wind or ACE 
near Earth), and it was found that the observations of- 
ten differed greatly, even for small longitudinal separa- 
tion of a few degrees (see the e xtensive summary o f these 
recent and historic events by IKilpua et al.l 1201 lh . This 
might be expected for MFRs that have a high inclination 
of their axis to the ecliptic plane, because one space- 
craft passes through the center of the flux rope and the 
other one crosses along its e dge ([Mostl et al.ll2009aH b1). 
ICrooker fe Intriligatorl ([ 19961 ) showed that it is also pos- 
sible for a high inclination flux rope to have a strongly 
elongated cross-section. 

For flux ropes with a low inclination to the ecliptic, 
one can probe the distribution of the magnetic field along 
the tube, and significant differences of the axial orienta- 
tions and magnetic fluxes h ave been found at a space- 
craft separation of 20-40° ([Farrugia et al.l 1201 ID . pre- 
sumably caused_bj^he_interaction with a high-speed 
stream. IKilpua et al.l ([201 If ) summarized the available 
multi-point observations of ICMEs by noting that "the 
characteristics of ICMEs and the structure of the so- 
lar wind they are embedded in varied significantly from 
event to event" . 

Thus, the common picture of flux ropes inside ICMEs 
as magnetic loops extending from the Sun into interplan- 



etary space likely needs to be revised to include distor- 
tions of the axis and the cross-section. In fact, the mag- 
netic field geometry itself is also uncertain; while it is gen- 
erally thought to consist of helical field lines, other config- 
urations can lead to similar appearances of the magneti c 
field in situ ([Jacobs et al.ll2009t lAhHaddad et al.ll2Qlll ). 
so this matter is also open to further inquiry. 

Throughou t this paper, we stick to definitions given 
bv IRouillardl (l20lll) . We define an "ICME" as the full 
interval of solar wind signatures which has changed as 
compared to the background wind due to the presence 
of the eruption. This includes the shock, the sheath re- 
gion of piled- up and compressed plasma behind the shock 
and the magnetic ejecta, often driving the shock at 1 AU. 
Subclasses of ejecta are MFRs, which are defined by an 
interval of smoothly rotating magnetic field in combi- 
nation with enhanced total field strength. If the inter- 
val additionally fulfills low proton temperature they ar e 
called "magne t ic clo uds" (MCs; iBurlaga et "all 119811) . 
iLepping et al.l ([2006D also defined "magnetic cloud like" 
(MCL) structures which have the general characteristics 
of MCs, but their field rotation is not clear enough to put 
them into one of the e ight c atego ries for MCs defined b y 
Bothmer fe Schwennl (119981) and iMulligan et~aH (l998). 



Mostl et al.l (|2010l) and iWood et alj (|201lD have found 



a case where the MCL signature was definitely caused 
by the spacecraft crossing a leg o f an ICME seen in he- 
liosph eric images, as assumed bv iMarubashi fc Leppingl 

(pool . 

All the above ideas are made even more complicated by 
the fact that coronal mass ejections (CMEs) often tend to 
erupt in rather quick succession, a feature which was orig- 
inally obse rved for flares and thus called "sympathe tic" 
flaring (see lSchriiver fc TithlfflTTl : IToTok et al.1120111) . In 
such a case, one eruption is causally related to the other. 
In the interplanetary medium, given these favorable ini- 
tial conditions, ICMEs might run into one another and 
interact. Especially around solar maximum, with the 
next one expected for 2013, this could be a rather com- 
mon situation (e.g., recent periods of strong successive 
activity include 2011 February and 2012 March). If as 
a result of the merging process the identification of the 
individual ICMEs is impossible, the st ructure is called a 
"complex ejecta" ([Burlaga et al.|[2002T) . However, it may 
also happen that coalescing events retain their individ- 
ual structure, but the mutual interaction involving en- 
ergy and momentum transfer between the ICMEs may 
alter their field and plasma structu r e in distinct ways 
(e.g.. iFarrugia fe Berdichevskvl 12001: iLugaz et al.l 120091: 
ILiu et al.H2012T) ~ In general, solar wind streams which 
show two speed maxima (either CME-CIR or CME 
CME) are also called "compound streams" and they 
are associated freque ntly with large geomagnetic storms 
([Burlaga et al.l fl987l ) . Thus, these streams are of great 
interest to space weather research in general, simply be- 
cause the disturbed solar wind interval is both longer and 
compression may lead to enhanced and m ore geo-effective 
southward pointing magnetic fields (e.g., iDal Lago et al.l 
120061 : iZhang et al.ll2007E iRouillard et al.ll2010ft 

A series of CMEs took place on 2010 August 1, asso- 
ciated with several flares, filament eruptions, and coro- 
nal dimmings, and was one of the first periods of strong 
activity in the rising phase of the current solar cycle 
24. In a very fortunate way, many spacecraft equipped 
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with suitable instruments and far separated in heliolon- 
gitude were positioned at different heliospheric distances 
to study this event chain from Sun to Earth with both 
imaging and in situ observations. To get the most out 
of this unprecedented dataset, an international team has 
formed to discuss the events at workshops in Abingdon 
(UK), Graz (Austria) and Aberystwyth (UK), all held 
in the first half of 2011. Our paper belongs to a se- 
ries of studies that form the output of these workshops 
in an attempt to characterize these events thoroughly, 
and to lay the foundation for future studies on these 
events so to obtain more new insights into the physics 
of solar eruptions. The solar coronal signatures revealed 
by the Atmospheric Imaging Assembly and Helioseis- 
mic and Magnetic Imager in s trume n ts on SDO were 
studie d in detail by ILiu et all (|2010ft . iSchriiver fc Title! 
(|2011ft and iLi et al.1 (|2011| ). The associated CMEs and 
their connection to the ICMEs o bserved at Earth us- 
ing the Heliospheric Imager (HI; lEvles et al.ll2009f) on- 
board STEREO- Ah e ad ha s been made bv lHarrison et al.l 
(|2012ft . iWebb et al.l (|2012D linked th e flows seen with the 
Solar Mass Ejection Imager (SMEI. Wivles et al.l[2003h to 
in situ observations of MFRs. Our paper is meant to 
close the circle, and to analyze all the available in situ 
observations in detail, as part of the most complete data 
set on a sequence of solar eruptions we are aware of, ob- 
served from the Sun end-to-end to 1 AU. 

The available in situ observatories were spread around 
the inner heliosphere at the beginning of 2010 August in 
an almost laboratory-like configuration, as seen in Fig- 
ure [Ha). Venus Express (VEX), in a polar orbit around 
Venus at the time of the event, was situated between 
STEREO-Behind and the Earth, and the MESSENGER 
probe was located between STEREO-Ahead and the 
Earth. Thus this spacecraft configuration covered 150° in 
heliospheric longitude. The event also happened before 
MESSENGER 's insertion to an orbit around Mercury in 
2011 March. The observations and analysis we carry out 
sets strong constraints on the ICME's propagation direc- 
tions and orientations as well as their arrival times, which 
can serve as a benchmark test for space weather pre- 
diction efforts usin g numerical simulations (carried ou t 
for these events bv lOdstrcil et alj|2012t IWu et al.|[20Tlh . 
They are also of great value in disentangling the signa- 
tures of the first multipl e ICME fronts obse rved end- 
to-end by STEREO/HI (jHarrison et all 120 121 ) . and put 
constraints on the kinematics of a CME-CME interac- 
tion event (ITemmer et all2012tlMartinez- Oliveros et al.l 
I2012t ILiu et alj|2012f >. In this sense, we think that this 
series of events, observed at extremely high resolution 
on the Sun by SDO, followed to 1 AU with STEREO/HI 
and SMEI, and subsequently observed in the solar wind 
at multiple locations may serve as a standard event for 
many future investigations on multiple and interacting 
ICMEs. 

2. SPACECRAFT POSITIONS AND OVERVIEW OF 
OBSERVATIONS 



10°, with other direction-finding techniq ues (see also 
ITemmer et"all l2012t lHarrison et ail 12012ft . From Fig- 
ure [Ha) ! where these directions are indicated as solid 
red arrows (a shorter arrow indicates a later eruption), 
it is seen that they are consistent with the longitudes of 
the in situ spacecraft, so they can be expected to leave 
ICME signatures at all of them, aside from STEREO- A, 
if they propagate more or less radially. Another CME on 
2010 July 30 (the dashed arrow) was directed mainly to- 
wards STEREO-B, located at E71. This previous CME 
is of importance because it influenced the prop agation of 
the A ugust 1 CMEs in the eastern direction (|Liu et al.l 
1201% 

FigureQJa) shows the spacecraft positions on 2010 Au- 
gust 3 00:00 UT, summarized in Table [1] as well as the 
initial CME directions and speeds, and structures ob- 
served in situ at each spacecraft. The largest separation 
in heliospheric longitude between STEREO-Behind ( "B" 
in the figure, at a radial distance from the Sun of 
r B = 1.060 AU, and most eastward) and STEREO- 
Ahead ("A", rA — 0.961 AU, most westward) was 
149.86° in Heliocentric Earth Ecliptic (HEE) coordi- 
nates. To simplify the analysis we further neglect the 
small separation angles in HEE latitude between the 
spacecraft. The MESSENGER (r M = 0.38 AU) and 
VEX [rv = 0.727 AU) probes were in extremely favor- 
able positions, each roughly halfway between the respec- 
tive STEREO spacecraft and Wind, which was located at 
the LI point of the Sun-Earth system (rw = 1.006 AU) 
at the time. Note that MESSENGER in its close orbit 
to the Sun moves much faster than the other spacecraft. 
From August 1 to 4 r^j changes from 0.40 to 0.37 AU 
and the longitude from 42.9° to 52.6° . 

Figure [TJb) shows from top to bottom the to- 
tal magnetic field strength |B| at the locati ons of 
STER EO-B (lAcuna et al.l 120081: ILuhmann etHI [2008V 
VEX (iZhang et . all 12006ft. Wind (iLepping et all 1*1995). 



IQdstrcil et "all (|2012ft identified four significant CMEs 
on 2010 August 1 with initial ( < 15 i? ) direc- 
tions ranging in longitude from 25° east to 8° west 
with respect to the Sun-Earth line. We further quote 
their results from an elliptical forward modeling tech- 
nique, which are in general consistent, mostly to within 



MESSENGER (lAnderson et all 120071). an d STEREO-A 
(|Acuna et aT]|2008l : ILuhmann et al.H 2008). ordered from 
top to bottom by heliospheric longitude, from east to 
west. Henceforth we designate by Mi (i=l, 2, 3) large- 
scale magnetic flux ropes and by Si (i= 1, 2) the shocks 
(to be discussed in detail in the later sections). 

For the moment we do not imply any associations be- 
tween the structures observed at each spacecraft. This 
plot immediately shows that the |B| profile observed by 
VEX is qualitatively similar to |B| at STEREO-B, at 
18° separation in longitude, while it looks quite different 
at Wind and MESSENGER. At the latter two space- 
craft, a single shock is present, in contrast to two clear 
shocks at STEREO-B and VEX. The featureless pro- 
file at STEREO-A does not show signs of any shocks or 
ejecta. 

2.1. Timing Considerations 

It is tempting to think that the shocks SI and S2 
in Figure [TJb) are signatures of two different, extended 
shock fronts. By this we mean that there actually exist 
only two large-scale fronts, which the spacecraft sam- 
ple at different points in space and time. We can jus- 
tify this by estimating the arrival times of the shocks if 
the spacecraft were all positioned at the same heliocen- 
tric distance, in our case the heliocentric distance of the 
Wind spacecraft (rw= 1-006 AU), which we further use 
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(a) Spacecraft positions, CME directions and ICME structures (b) 
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Fig. 1.— (a): Positions of spacecraft and inner planets on 2010 August 3 00:00 UT in the ecliptic plane (B= STEREO-B, A=STEREO-A). 
The principal directions of one CME on July 30 (das hed) and four CMEs on August 1 (so lid, numbered 1-4), are shown by orange arrows, 
with a shorter arrow indicating a later eruption (Harrison ct al. 2012; Odstrcil ct al. 2012). We also indicated the initial de-projected CME 
speeds (Odstrcil ct al. 2012) for distances < 15 Rq. For each spacecraft, the progression of shocks (S) and flux ropes (M) is indicated, 
(b): The solar wind total magnetic field observed at five locations in the inner heliosphere (< 1 AU) close to the ecliptic plane. From east 
(STEREO-B) to west (STEREO-A), the covered heliosphcric longitude is 149.8°, and ICME signatures are observed from STEREO-B to 
MESSENGER, corresponding to 120.5°. From top to bottom (or east to west): the total magnetic field at STEREO-B, Venus Express 
(in orbit around Venus), Wind (at the LI point), MESSENGER, and STEREO-A. Vertical solid lines indicate shock arrival times, and 
dashed red vertical lines estimated shock arrival times extrapolated to the heliospheric distance of Wind (rw = 1.006 AU). Shaded intervals 
indicate well-defined magnetic flux ropes. 



TABLE 1 
Spacecraft Positions 



Spacecraft 


Abbreviation 


r (AU) 


r (Rq) 


HEE Long.(°) 


HEE Lat. (°) 


Bmax 


R' 

max 


STEREO-Behind 


B 


1.060 


228.1 


-71.29 


-0.29 


33.2 


36.2 


Venus Express 


V 


0.727 


156.2 


-53.18 


-0.04 


58.5 


34.4 


Wind 


W 


1.006 


216.4 


0.16 


-0.03 


19.2 


19.2 


MESSENGER 


Al 


0.38 


81.7 


49.2 


-5.3 


33.4 


6.8 


STEREO-Ahead 


A 


0.961 


206.7 


78.57 


0.01 


10.1 


9.4 



Note. — Quoted for 2010 August 3 00:00 UT. The longitude and latitude are given in Heliocentric Earth Ecliptic (HEE) coordinates, 
with the Earth at 0° in both latitude and longitude and the longitude being positive to the west of the Earth. Bmax is the observed 
maximum magnetic field strength between 2010 August 1-5, and B' max the scaled field strength to the heliocentric distance of the Wind 
spacecraft (see the text). 



as a reference point. 

To this end, we plot the extrapolated arrival times at 
rw of both shocks as red vertical dashed lines in Fig- 
ureQJb), by using constant speeds of Vsib — 460 km s _1 
for SI, and V S2 b = 600 km s" 1 for S2 at STEREO-B, 
which are the shock normal speeds at STEREO-B (de- 
rived in Section 3). For VEX there is no high resolu- 
tion plasma data available, so we assume this speed is 
Vs2V = 550 km s -1 , a mean value between the results 
for Wind and STEREO-B. While the initi al speed of the 
most eastward CME2 is « 1200 km s -1 . iTemmer et al~l 
(|2012f) showed that there was a strong deceleration due 
to interaction with CME1 around 30-50 Rq, and that 
CME2 attained a constant speed well before the helio- 
centric distance of Venus. It is also not expected that 



the IC ME speed chang es much between 0.72 AU and 
1 AU (jJian et al.l I^OPS) . Thus, assuming this speed to 
be constant from the Venus distance to about 1 AU 
and beyond is reasonable. For MESSENGER we use 
Vs2M = 500 km s _1 , similar to the value observed for 
Wind. Again, for MESSENGER there is no high-time 
resolution plasma data. If we assume some deceleration 
up to 1 AU, this choice is consistent with the ~ 600 
km s _1 initial speeds fo r both CMEs head ing to the 
west of Earth (CME3/4, fOdstrcil et aTll2012[ ). and sim- 
ilar propagation speeds for these CMEs derived from 
STEREO/HI observations bylflarrisoneLalJ ((2012J). We 
also assume that the shock surface moves radially away 
from the Sun. 
Figure QJb) shows a clear progression from east to 
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west for the time-shifted arrivals (see the red dashed 
lines): S2 arrived first at STEREO-B, then at VEX, 
followed by Wind (all three within 12 hr), and last at 
MESSENGER, and it is undetected by STEREO- A. The 
extrapolated arrival times at STEREO-B (vertical red 
dashed lines) in Figure [1] are earlier than the observed 
ones because STEREO-B is at a larger heliocentric dis- 
tance than Wind. The close timing compels us to be- 
lieve that we are observing the same shock surface at 
the different locations, and the late MESSENGER ob- 
servations and the lower shock normal speed at Wind 
are expected because for a convex surface, such as an ex- 
panding circular front attached to the Sun, the side parts 
(or "flank") will move slower than the fro nt part (also 
called "apex" or "nose" of the ICME; cf. IMostl et all 
120111: IMostl k Daviesl [2012h . We thus deduce qualita- 
tively that the shock surface for S2 has a convex shape, 
with a smaller radial distance of the shock front corre- 
lated with positive longitude values for a given point in 
time. This means that the shock front is in this case not 
spherical with respect to the Sun. 

There are two effects which we think contribute to this 
shape. First, the previo us ICME on Ju ly 30, which is as- 
sociated with shock SI (Liu ct al. 2012|), created a region 
of depleted density in the solar wind, so the eastern part 
of the ICMEs originating on August 1 experienced less 
drag, which led to less deceleration. Secondly, one has to 
keep in mind that the initial speed of CME2, with an ini- 
tial direction between E25 and E41, was » 1200 km s _1 
and thus about double than for the other CMEs (all 
» 600 km s- 1 ). 

The shock SI is constrained to the east of Earth and 
is observed first at STEREO-B, then at VEX, but nei- 
ther at Wind nor MESSENGER so we cannot derive an 
approximate global shape like for S2. 

2.2. Shape of the Shock Fronts 

We have also deduced the shape of the shock front 
in a more quantitative way. In many situations, this 
shape is unknown or accessible only through numerical 
simulations. However, this shape is of great interest for 
predicting the arrival tim e of ICMEs as a functi on of he- 
liospheric longitude (e.g... IMostl k Da vics 2012|), and we 
are especially interested in its appearance for a complex 
situation involving interacting ICMEs. 

In Figure [2ja) we plotted the estimated heliocentric 
distance of the shock (e.g., rs2B for the shock S2 at 
STEREO-B) along radials from the Sun to each space- 
craft for the time the shock is observed at Wind (tsiw 
= 2010 August 3 17:05 UT) using 

rS2B = r B +Vs2B(ts2W - *S2b)) (1) 

with Vs2B and ts2B being the S2 shock normal speed 
and its arrival time, respectively at STEREO-B. Again, 
we assume radial propagation of the shock. The result 
is rs2B = 1-2349 AU, which is considerably larger than 
r w = 1.006 AU. Both distances would be similar if the 
shock had the shape of a sphere centered on the Sun. 
Similarly, estimating these distances for VEX and MES- 
SENGER by using similar speeds as for the time-shifts 
(550 km s- 1 for VEX and 500 km s" 1 for MESSEN- 
GER), we get r S2 v = 1-1188 AU and r S2M = 0.7412 AU. 
We may thus plot the shock front S2 as the red solid line 



connecting the dots (cf. iBurlaga et al.lH98ll ) at the snap- 
shot of the shock arrival time at Wind in Figure [2f a), 
showing clearly the non-spherical shape. Between VEX 
and Wind as well as Wind and MESSENGER we plotted 
another point, estimated by eye for a smoother form of 
the shock surface. 

A similar analysis was done for shock SI, which seems 
to have a shape closer to spherical than S2, but with 
only two data points it is difficult to discuss the global 
structure of SI so we will further concentrate on S2. 

As we will discuss in a later section, we note that if 
S2 is indeed one large shock surface spanning > 120° in 
heliolongitude, it likely is driven by different structures 
at different locations along the front. 

2.3. Scaling the Magnetic Field Strength 

Additionally, to fully characterize the observed com- 
plex flows in a simple fashion, we can deduce how the 
maximum observed magnetic field strength scales with 
heliospheric longitude, and again, we need to roughly 
compensate for the different heliospheric distances r of 
the spacecraft. As a proxy, we use an empirical scalin g 
relation for the inner heliosphere bv lLeitner et al.l ([2007) , 
who found that the maximum magnetic field strength Bq 
inside MCs in the inner heliosphere (< 1 AU) depends 
on radial distance as 

B (r) = lS.lr- 1 - 64 , (2) 

with Bq in (nT) and r in (AU). For example, the peak 
field strength was observed by VEX at the front of the 
flux rope M2 as B max y = 58.5 nT at a radial distance 
of ry . To estimate the maximum B' maxV scaled to the 
radial distance of rw, we simply eliminate the factor 18.1 
to obtain 

B' maxV = B maxV (r w /r v )~ 1SA , (3) 

yielding B' maxV = 34.4 nT. At STEREO-B and MES- 
SENGER, the peak field (see Table [TJ is also detected 
shortly after S2, while it is more delayed to S2 at Wind, 
our reference location. Taken together, the maximum 
field strengths estimated at r\y compare as: 36.2 nT 
(STEREO-B), 34.4 nT ( VEX), 19.2 nT {Wind, not 
scaled), and 6.8 nT (MESSENGER). For comparison, 
the result is 9.4 nT for STEREO-A, which did not ob- 
serve the structure. Thus, the high peak field at MES- 
SENGER is scaled down to normal solar wind levels at 
1 AU. This also means that the shock, if it had an initial 
component towards STEREO-A, is likely to have decayed 
when it reached 1 AU at its western edge. 

Because one gets easily lost in these huge amounts of 
data, it is useful to keep the basic picture which emerged 
from this section in mind as a background for the follow- 
ing, more detailed analysis. We now move on to calcu- 
late parameters of the shocks and flux ropes to obtain 
an even more detailed global picture using the in situ 
observations. In addition to parameters discussed so far, 
these may be useful for more detailed comparisons to 
numerical simulations. 

3. SHOCK ANALYSIS 

The aim of this section is to analyze the observed 
shocks in detail and to relate these findings to the global 
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(a) Global shock surface from in situ observations (2010 August 3 17:05) 
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Fig. 2. — (a) Global configuration of two shock fronts SI and S2, estimated at the time of the shock observation at Wind (2010 August 
3 17:05 UT). The approximate shape of the shocks (red solid lines) is shown by estimating the heliocentric distance of the shocks at each 
location for the given point in time. For SB, shock normals are drawn as arrows in the color consistent with the spacecraft position, the 
configuration of the shock is indicated as oblique or quasi-perpendicular. The nominal Parker spiral field is drawn connecting the Sun 
to each spacecraft for the upstream conditions of solar wind speed 450 km s _1 at STEREO-B and VEX and 400 km s — 1 at Wind and 
MESSENGER. Note that VEX and STEREO-B are approximately connected by the nominal Parker spiral field, (b) The insert shows an 
illustration of the four-spacecraft method to determine the shock normal orientation near Earth, using Wind, ACE and ARTEMIS P1/P2. 
A planar shock is measured at four points in space for which the distances to one reference spacecraft have to be known. 



TABLE 2 

Shock Parameters at Each Spacecraft. 



Spacecraft 


Shock 


Long.(°) 


Arrival (UT) 


Normal n s (xyz) 


8s(°) 


4>s{°) 


V s (km s- 1 ) 


<z(°) 


M A 


Bd/B u 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


STEREO-B 


SI 


-71.3 


Aug 2 15:30 


(-0.998. 0.058. -0.004) 





177 


460 


69 


1.9 


2.3 


VEX 


SI 


-53.2 


Aug 1 14:41 


(-0.761, 0.263, -0.594) 


-36 


160 




56 




2.3 


STEREO-B 


S2 


-71.3 


Aug 3 05:00 


(-0.981. 0.194. 0.025) 


1 


169 


600 


79 


3.1 


1.7 


VEX 


S2 


-53.2 


Aug 2 11:30 


(-0.981, -0.094, 0.182) 


10 


185 




54 




2.6 


Wind 


S2 


0.2 


Aug 3 17:05 


(-0.831. -0.259, 0.492) 


29 


197 


497 


56 


4.6 


3 


Wind, ACE, ARTEMIS P1/P2 


S2 


~ 


Aug 3 17:05 


(-0.907, -0.218, 0.361) 


21 


194 


511 


42 


5.0 


2.6 


MESSENGER 


S2 


49.2 


Aug 2 09:43 


(-0.794, 0.345, -0.501) 


-30 


157 




44 




1.4 



Note. — Column 1: observing spacecraft. Column 2: name of the shock. Column 3: longitude of the spacecraft in HEE coordinates. 
Column 4: shock arrival time. Column 5: shock normal in GSE coordinates for near-Earth spacecraft, and approximated GSE coordinates 
for the other spacecraft (i.e., RTN rotated by 180° in the RT plane for MESSENGER and STEREO-B and VSO coordinates for VEX). 
Column 6: latitude of the shock normal angle with respect to the ecliptic plane. Column 7: longitude of the shock normal measured from 
GSE +X (0°) to GSE +Y (90°). Column 8: shock normal speed. Column 9: angle between the shock normal and the upstream field (also 
known as Bsn)- Column 10: Alfven Mach number. Column 11: ratio between the downstream and upstream field strengths. 
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shape and structure of the shocks, giving us a rare oppor- 
tunity to d o this on a global sc ale with multiple space- 
craft (e.g., iBurlaga et al.l IT981). Parameters of inter- 
planetary shocks, such as the shock normal orientation, 
the configuration of the shock (perpendicular, oblique 
or parallel), and the Mach number, can be obtained by 
using the co-planarity theorem, which follows from the 
general Rankine-Hug oniot jump conditions of magneto- 
hydrodynamics (e.g. , iKivelson fc Rus sell 199H iBurlagal 
119951: iSchwartzl 119981 ). The theorem states that the up- 
stream (index u, or "unshocked" or "undisturbed" re- 
gion of the solar wind) and downstream (index d, or 
"disturbed" ) vectors of the interplanetary magnetic field 
(IMF) and the shock normal n s all lie in the same plane, 
and the normal can be obtained with 



(B u - Brf) x (B u x B d ) 
|(B U -B d ) x (B u xB d ) 



(4) 



For near-Earth spacecraft, we will express the result- 
ing shock normal in geocentric solar ecliptic coordinates 
with X pointing toward the Sun, Z toward ecliptic north 
and Y completing the right handed triad (thus point- 
ing towards solar east). Also, we will quote angles for 
the shock normal 9 S (latitude) and 4> s (longitude) which 
are defined as such that 9 S is measured with respect to 
the ecliptic plane (9 S = 0°, and ecliptic north = +90°) 
and (f> s is 0° along GSE X and 90° along GSE Y. For 
STEREO- A/B and MESSENGER, we will simply rotate 
the result, first given in radial-tangential- normal (RTN) 
coordinates, which are defined as R pointing away from 
the Sun and T towards solar west, by 180° in the solar 
equatorial plane, meaning we replace R with —R and T 
with — T, so the orientation of the coordinate system is 
similar to GSE. For VEX, Venus Solar Orbital (VSO) co- 
ordinates are used, which are defined similar to GSE but 
for the orbit of Venus. We are thus neglecting small dif- 
ferences in orbital inclination between the Earth, Venus, 
and the solar equatorial plane (see also Table [l}. Typical 
errors for the shock normal with the co-planarity method 
are in the order of 10°-20°, so the intrinsic errors are def- 
initely larger than the systematic error we impose by our 
simplified treatment of the coordinate systems. 

If solar wind bulk velocity vectors are available (only 
for the Wind spacecraft), one can use the conservation 
of mass flux across the shock surface to obtain the shock 
normal speed V s with 



Ys 



N d V d 



N d -N u 



(5) 



where N d , N u , Vd, and V u are the downstream and 
upstre am solar wind densities and velocities, respec- 
tively (jKivelson fc RusselH H995). In the case when no 
plasma bulk velocity vectors are available but its norm 
is known (for STEREO- A/B), we can still estimate the 
shock speed if we suppo se that upstr eam and down- 
stream speeds are radials (|Burlagalll995l ) and replace V u 
and Vd by their respective norms V u and V d in Equa- 
tion ©. For VEX and MESSENGER we will derive 
the shock normal only because no high resolution plasma 
data are available. 

The Alfvenic Mach number defines the strength of the 
shock and it is calculated as Ma = \V sr — V u \/Va, where 
Va is the Alfven speed in the upstream region and V sr the 



shock radial speed (V sr = V s /n sx ). The angle q (often 
called 9 Bn in the literature, but we would like to avoid 
this because in this paper 9 denotes latitudes for flux 
rope and shock orientations) , gives the angle between the 
upstream field B u and n s and defines the shock as quasi- 
parallel when q < 45° and quasi-perpendicular if q > 45°. 
Shocks with (/-values around 45° are called oblique. The 
parameter q d efines the mec hanisms of particle accelera- 
tion at shocks (Reames 1999), and sh ocks propagating in- 
side I CMEs have higher values of q (Richards on fc Cand 
120101 ). At 1 AU, for a radially propagating spherical 
shock, q is expected to be close to 45° because it is the 
angle of the Parker spiral field to the radial direction. 

3.1. Application to STEREO-B, VEX and 
MESSENGER 

We will now discuss the application to the STEREO - 
B shocks and summarize briefly the results at VEX and 
MESSENGER. Figure [3] shows the p lasma and magnetic 
field measure ments by PLASTIC ([Galvin et al.l I2008T) 
and IMPACT (|Luhmann et al.l l2008f) aboard STEREO- 
B. 

Shock SI : We used an interval August 2 05:08- 
05:53 UT. The results for this method are n s u, — 
(-0.998,0.058,-0.004) (in coordinates where R is re- 
placed by -R and T by -T), or 9 slB = 0° / <p s i B = 177°, 
and q\B = 69°. The Alfvenic Mach number was 1.9 and 
the ratio of the upstream to downstream field was 2.3, 
both values pointing to a moderate to weak shock. 

Shock S2: For interval August 3 04:40-05:20 UT 
we get n s2B = (0.981,-0.194,0.025), or 9 s2B = 1° / 
4>s2B — 169°, and q 2B = 79°. Here, Ma was 3.1 and 
the B d / B u =1.7, again in a moderate regime of shock 
strength. For a compilation of these results, see Table [2j 

Figure |4] shows the magne tic field measured by 
VEX/MAG (|Zhang et all 12006). We linearly interpo- 
lated over intervals when VEX was closer than two 
Venus radii to the planet to avoid plotting field strengths 
strongly exceeding those in the solar wind when the 
spacecraft flies through Venus' induced magnetosphere. 
The application of the co-planarity theorem shows that 
both shock normals were close to radial from the Sun. 
The shocks were of oblique type, and of moderate to 
weak shock strength. 

MESSENGER observed a single shock on August 2 
9:43 UT as seen in Figure [5j The derived parameters 
indicate a very weak shock strength with B d /B u = 1.4 
and the strongest deviation of the shock normal from 
the radial direction. The low strength at this western 
location could also indicate that if the shock S2 had an 
initial component directed towards STEREO-A, it could 
have already decayed when it reached this spacecraft, 
as noted previously when we estimated how the field 
strength scales with heliospheric distance. 

3.2. Near- Earth Space 

Wind is currently in orbit around the Lagrangian 
point LI and at the time of the event it was located 
at R w = (205, -64, -14) Earth radii (R E ) in GSE coor- 
dinat es. Magnetic field (M FI instrument: iLepping et al.l 
[1991 and plasma (SWE; lOgilvie et all [1995D obser^aT 
tions between 2010 August 3 and 4 by the Wind space- 
craft are shown in Figure [6l Because we are interested 
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STEREO-B IMPACT/ PLASTIC 2-4 August 2010 




_l I l l I I i i L_ 

12 18 00 06 12 18 00 06 12 18 



03 Aug 04 Aug 



Fig. 3. — From top to bottom: magnetic field components —Br (red), —Bt (green), and total field magnitude (black); magnetic field 
B z component; proton bulk speed; proton number density; proton temperature; and proton /3. Two shocks are indicated by vertical solid 
lines along SI/ 2 and magnetic flux ropes are indicated by shaded regions M 1/2/3. 
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mainly in the relation of the magnetic field vector to the 
same observations by the other widely separated space- 
craft VEX and STEREO-B, we plot here the magnetic 
field components in GSE coordinates. Later, when dis- 
cussing geoeffects, we will use the appropriate GSM co- 
ordinates. 

At 20f0 August 3 f7:05 UT, Wind detected a forward 
shock that is indicated in Figure [6] with a solid line and 
labelled S2. At the shock the solar wind speed increased 
abruptly from 425 km s^ 1 to about 505 km s _1 . Simul- 
taneously, the ion density increased from 5 to 12 cm~ 3 , 
and the magnetic field magnitude from ~3 nT to ~9 nT. 
Note that this shock shows, when zoomed in, a two- 
step profile, so the solar wind speed actually increases 
to ss 600 km s^ 1 later in the sheath region. Here, in 
addition to the co-planarity theorem we were able to use 
a four-spacecraft method which is based on the different 
shock arrival times and the locatio ns of the spacecraft 
(jRussell et al.lH983t lSchwartzlU998l ). The shock S2 ar- 
rived at A CE on August 3 1 6:55 UT, and at the two 
ARTEMIS (|Sibeck et al.|[20TTl former THEMIS) space- 
craft PI on 17:20 UT and P2 at 17:40 UT, while both 
were in transit from Earth to lunar orbit. 

Method (1): 

We used 3 s magnetic field measurements from 
Wind in the interval 16:50-17:20 UT. The co-planarity 
approach yielded a shock normal orientation n — 
(-0.831,-0.259,0.492) in GSE coordinates (latitude 
6 S = 28°; longitude <j) s = 197°). The shock speed 
was V s = 497 km s _1 , and the shock was oblique with 
q = 56°. From Va = 38 km s _1 being the Alfven speed 
in the upstream region and V u = 410 km s _1 the up- 
stream flow speed, we get an Alfven Mach number of 
about 4, and the ratio of the downstream to upstream 
field is Bd/B u = 3. Both values point to a moderately 
strong shock. 

Method (2): 

Figure[2][b) shows an illustration of the four-spacecraft 
method. It assumes local planarity of the shock and a 
constant propagation direction and speed. When the 
spacecraft positions i?; and the differences At a i (with 
a = 2,3,4) of the shock arrival times are known, the 
shock normal n s and speed V s can be obtained from: 

V s = v ■ n s , (6) 
(Rq, - Ri) • n s = v ■ At ai , (7) 

with v being the norm of the shock normal speed. The 
results for this method are n s = (-0.907, -0.218, 0.361), 
or orientation angle 9 S = 21°, (f> s = 194°, and q = 40°, 
V s = 511 km s _1 and Ma = 5.0. We thus find a 
good agreement between this and the single-spacecraft 
method, with differences in the shock normal of only 8° 
and in the speed of 14 km s _1 . This means that the 
shock normal near Earth is well determined, indicating 
that the shock surface is reasonably close to planar on a 
scale of sslOO Re- 

3.3. Implications for the Global Shock Shape and 
Structure 

We now discuss the individual shock parameters, sum- 
marized in Table [21 in view of the question: may we 
extrapolate features of the shock surface from a local to 
the global scale? 



It was found in previous studies that even on a length 
scale of o nly a few Re, sho ck normals can be inconsis- 
tent (e.g., Sz abo et "all 120011 ). pointing to non-planarity 
or corrugation of the shock surface on very small scales 
compared to the global extent. We can observe the shock 
S2 at a longitudinal spacecraft separation of 120°, which 
corresponds to a length of l s ~ 2.1 AU at the heliocen- 
tric distance of 1 AU (neglecting the non-spherical shape 
of the shock). The maximum separation of ncar-Earth 
spacecraft is of the order of a few hundred Re, which is 
roughly a factor 100 less than l s . It thus becomes clear 
that for a corrugated shock surface on scales much less 
than l s , it would be impossible to extrapolate the global 
shape of the shock from a sin gle-spacecraft observation. 
However, iSzabo et al.l (|2001[) also note that fast ejecta 
impose a more pronounced planarity to the shocks which 
they are driving as compared to slower ones. A prime ex- 
ample for such a "well-behaving" s hock near Eart h drive n 
by a fast ejecta is discussed by iFoullon et al.l (|2007j) . 
It may thus be expected that, for a fast compound 
stream like the structures following S2, with speeds up 
to 800 km s" 1 at STEREO-B and 600 km s" 1 at Wind, 
it could indeed be possible to extrapolate the principal 
direction of the shock front from the single-point shock 
normals and see if it is consistent with the "real" situa- 
tion as determined by the multi-spacecraft observations. 

To this end, we may take a look again at Figure [2f a). 
In particular, note the longitude of the shock normals, 
which flare away from the radial direction, towards the 
east at STEREO-B, slightly west at VEX and more pro- 
nounced to the west at Wind. This is intuitively ex- 
pected for a non-corrugated shock surface which is con- 
vex with respect to the Sun. Only MESSENGER does 
not follow this pattern, which could be a hint that the 
shock surface at the outer flank is not driven strongly 
enough. An observer looking at near-Earth spacecraft 
would thus have correctly deduced, by looking at the lon- 
gitude of the shock normal, that the apex of this shock 
is to the east of Earth. This means that the principal di- 
rection of the shock can be correctly obtained in this case 
from single-point (using co-planarity) or multi-point (the 
four-spacecraft time-of-flight method) shock analysis of 
a shock driven by a fast ejecta in near-Earth space. 

The configuration of the shock S2 seems to be altered 
by the presence of the previous ICME at the location 
of STEREO-B, where the shock strength is l ower, due 
to th e higher field strength and lower density (|Liu et al.l 
120121 ) . Otherwise, the strength of the shock is moderately 
strong for the full extent. The shock S2 is also every- 
where oblique except at STEREO-B, where it is quasi- 
perpendicular, which is consistent with expectations. In 
general, shocks are more likely to be quasi-perpendicular 
while propagating inside an ICME (Richardson fc Cand 
2010). At Wind and VEX, the shock has just reached 
the back of a previous ejecta and seems not yet strongly 
altered in its configuration. From both methods we used 
in near-Earth space, the angle q is still close to 45° , as ex- 
pected for a spherical shock propagating through a nom- 
inal Parker spiral field (cf. Figure [2fa)). In contrast to 
the shock normal, the shock configuration is a function 
of the upstream field, which complicates the situation. 
Nevertheless, for this event it looks like the configura- 
tion of the shock as a function of longitude behaves as it 
is expected from our general understanding of interplan- 
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Fig. 4. — From top to bottom: total magnetic field and field components observed by VEX/MAG, in Venus Solar Orbital (VSO) 
coordinates (similar to GSE, but for the orbit of Venus). The bottom panel shows the radial distance of the spacecraft to the center of 
Venus. We linearly interpolated the magnetic field over times VEX approached Venus closer to two Venus radii to see only the solar wind 
magnetic field. 
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etary shocks. 

In summary, the global shape of the shock S2 is clearly 
non-spherical, though convex with respect to the Sun as 
expected, and the longitude of the shock normals n s , fol- 
low a pattern consistent with this shape at three out four 
locations. The following relation is valid: if the longitude 
of n s points east (west) of the radial direction, the shock 
apex is to the west (east) of the spacecraft. Only at 
the most western location at MESSENGER, which rep- 
resents its flank, the opposite is seen. This also places 
the apex of the shock between STEREO-B and VEX. 
The configuration of the shock as a function of helio- 
spheric longitude is complicated by the presence of pre- 
vious ICMEs, but follows a pattern expected for obser- 
vations at 1 AU. 

4. FLUX ROPE ANALYSIS 

In this section we use different modeling techniques 
to extract more information on the observed large-scale 
magnetic flux ropes (MFRs), i.e., those parts of the 
ICME intervals which exhibit a smoothly rotating field 
over a large angle, in addition to higher-than-average 
magnetic field. These are the structures we have de- 
fined as Ml/2/3 at the spacecraft STEREO-B, VEX 
and Wind, and Ml at MESSENGER. Our intention is 
to find matches between the model results to find out 
if the structures observed at a single spacecraft may be 
connected to each other and how the resulting global 
structures may look like. We use the same conventions 
for coordinate systems and orientation angles as for the 
shocks in the previous section. 

4.1. Flux Ropes Observed by Wind 

An inspection of Figure HI showing the magnetic field 
components in GSE coordinates and proton bulk param- 
eters, indicates that one MFR (labelled Ml) may be 
identified before the arrival of shock S2, and two dis- 
tinct MFRs can identified at Wind after the shock S2 
and its following sheath structure, which is the region 
of "shocked" solar wind with high temperature, density 
and proton j3 p . MFRs M2 and MS are immersed in a 
declining flow speed (V p ) pro file with two maxim a which 
is typical for complex ejecta (Burlaga ct al. 20Q2), and it 
also satisfies the definition of a compound stream. Thus, 
it is immediately clear that this flow is composed ofsey- 
eral structures, and, as was shown in iLiu et al.l {2012), 
these structures can still be treated as individual entities, 
also at the other spacecraft. Thus this stre am is not clas- 
sified as a "complex eject a" in the sense of lBurlaga et al.l 
(|2002f) . ILiu et al.l (|2012l ) suggested that structure Ml, 
whose back part seems to be "caught" by the shock just 
before arrival at 1 AU, could be the flank of a CME orig- 
inating on the Sun on July 30, so all ICME signatures 
except the proton temperature are not well pronounced, 
and thus we will not discuss this structure further here. 

The region M2 started on August 4 04:34 UT and it ex- 
tended almost seamlessly until the front boundary of MS 
arrived on August 4 10:07 UT, and MS extended further 
to August 4 23:18 UT. Both MFRs were associated with 
a low variance magnetic field and a smooth rotation of 
the magnetic field direction, and the proton temperature 
in MS is clearly low enough to classif y this structure as a 
magnetic cloud (|Burlaga et al.|[l98ll) . Also, (3 P showed a 



clear depression and the solar wind speed had a declining 
profile, the latter indicative of ongoing radial expansion 
when it passed Wind. These are typical signatures of 
MCs at 1 AU, thou gh its size of ~ 0.15 AU is not par- 
ticularly large (e.g.. Uian et alJ l2~006'). 

In contrast, the structure M2 shows stronger maxi- 
mum fields (19.4 nT compared to 14.1 nT for MS), higher 
proton temperatures, an enhancement in density denoted 
by Dl, a smaller size (~ 0.1 AU) and no radial expan- 
sion, which are all indicat o rs of i nteraction with MS after 
iFarrugia fc Berdichevskvl ((2004). It could also be that 
the shock S2 passed through M2 before 1 AU, and led 
to additional compression and heating, which is consis- 
tent with what is seen in heliospheric images (|Liu et al.l 
2012). At the leading edge of MS, there is another density 
increase (D2), which is stronger and of shorter duration. 
We do not think that it is the cause of compression, but 
rather that it is possible filament material. We plan to 
discuss this observation in detail in a subsequent publi- 
cation. 

More specifically, the average speed of the expansion 
is defined as a half of the difference between the speeds 
at the leading and trailing edges of the ICME. For MS, 
the measured speeds were 585 km s _1 and 504 km s -1 , 
respectively, yielding an expansion speed of 41 km s _1 . 
This is « 30 km s" 1 less th an typical for events observed 
at 1 AU (|Jian et al.l I2006D . On the contrary, M2 was 
associated with a positive speed gradient of 8 km s _1 
from the trailing to leading boundary, thus it was even 
contracting slightly as it passed 1 AU, which could also 
explain its smaller size. 

To estimate the local MC axis orientations and other 
parameters we used th e method of Grad-Shaf r anov (GS) 
recons t ruction (e.g.. IHu k. Sonnerud [2 002: Li u et al.l 
I2008at lM5stl et alj l2009alld : llsavnin et al.l 120111) . This 
method assumes time- independence and invariance along 
one direction which is equated with the axis of the MC. 
The first assumption is approximately fulfilled for M2 as 
shown above and also for MS since the radial expansion 
is in the order of 10% of the average speed for MS. By 
numerical integration, a map of the MC magnetic field 
lines in a plane perpendicular to the axis can be cre- 
ated without assuming a particular shape for the cross- 
section. The technique could be successfully applied to 
both MFR intervals. The resultant plots are shown in 
Figure [7] and the calculated values are again summarized 
in Table O 

Figures |7](a) and (b) show the transverse pressure 
P t = /2/i + p as a function of the vector potential 
A for each flux rope. Note that the plasma pressure is 
calculated from the protons only. For the assumed sym- 
metry, the vector potential may be expressed in the form 
A(x,y)e z , where z is the axis direction and B is then 
expressed as B = VA x e z + B z e z . The local MC axis 
orientation is found by constraining P t to be a single- 
valued function of A, i.e., having the same values along 
the in- and outgoing paths. To see how well this function 
can be fitted with a polynomial (of order 2, black solid 
line) with exponential tails (black dashed lines) which is 
used for the numerical i ntegration, a fitting residue Rf is 
defined (|Hu et al.ll200l . For M2, the value of R f = 0.06 
indicates a very good fit. For MS, this function has to 
be quite strongly extrapolated, and the fitting residue 
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Fig. 5. — Total magnetic field and magnetic field components observed by MESSENGER. For easier comparison to the components 
observed by the other spacecraft, component R is replaced by — R (pointing radially to the Sun) and component T by — T (pointing 
towards solar east). The solid vertical line marks the shock S2, and the shaded interval a magnetic flux rope (Ml). 
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(Rf — 0.16) is relatively high, which indicates that the 
method is not working well and that the magnetic field 
map and the results in Table [3] for this structure should 
be taken with care. This is also clear from looking at 
the magnetic field components, where it is not straight- 
forward to discern a clear full rotation of the field. Fig- 
ures [7f c) and (d) show the magnetic field cross-sections 
as contour lines. The shape of a M2 appears to be rather 
deformed, possibly as a result of the interaction with MS. 
The shape of the resulting cross-section of MS is remi- 
niscent of a glancing encounter; the strongest axial field 
is at the left border of the integration domain. 

In summary, the Wind observations of MS are likely 
to reflect a glancing encounter with a large M FR (e.g., 
iMarubashi fc Leppind [2007t iMostl et all 120101) . In con- 
trast, M2 shows a clear magnetic field rotation on 
roughly half the scale size (0.07 AU), which could be 
caused by a more central crossing of a smaller ICME 
and / or compression by MS. 

4.2. Flux Ropes at STEREO-B. VEX and 
MESSENGER 

We now briefly discuss large-scale flux ropes observa- 
tions at the other spacecraft. Modeling was not as suc- 
cessful as in near-Earth space, so we limit our discussion 
here. We also skip discussion of the flux rope Ml at 
STEREO-B and VEX, which follo ws shock SI and orig- 
inates from an eruption on July 30 (|Liu et al.ll2012j) . The 
main results are again summarized in Table 3. 

By looking at the magnetic field components in Figures 
[3] and [¥J we find very similar rotations for the structure 
MS, and at the relatively small separation of 18° between 
the two spacecraft it is quite clear that the same struc- 
ture is detected at two different sites. At STEREO-B it 
is a north-east-south (NES, right handed) flux rope with 
marginal rotation, classifying this as an MCL structure 
with an axial field pointing eastward and an axis close 
to the ecliptic plane. At STEREO-B, the speed profile is 
also roughly declining, as would be expected for an ex- 
panding flux rope. At VEX, one sees only the east-south 
(ES) part. Probably, this is a glancing encounter of an 
ICME flux rope which went to the north or south of the 
ecliptic. Note that non-force-free modeling gives consis- 
tent results in the orientation, yielding at both locations 
a low inclination to the ecliptic plane with eastward axial 
field (|Webb et alJl20Tl . 

At first glance, M2 at STEREO-B and VEX can be 
seen as a distinct MFR with a well-pronounced field ro- 
tation. It can also be successfully modeled as a small 
left-handed flux rope with moderate inclination pointing 
southward of the ecliptic. However, at both spacecraft 
there is no sheath region, so the rope follows the shock 
immediately, and thus we think it is much more likely 
that this is just the shocked part of the Ml structure 
and not a separate ICME. The scenario of the shock S2 
reaching the back of the CME on July 30 (connected to 
the structures SI /Ml) just before 1 AU is a lso consistent 
with heliospheric imaging (jLiu et al.ll20"H . 

At MESSENGER, we could identify only one short in- 
terval where the field is relatively smooth and rotating. 
This follows the shock S2 and is an NWS rope, with 
the axis being close to the ecliptic plane and the axial 
field pointing in the +T (RTN coordinates) or solar west 
direction. 



4.3. Plotting the Global Flux Rope Structure 

In Figure the local orientations of those structures 
for which we could find an axial orientation as given in 
Table [3] are visualized by different viewing angles. Fig- 
ure^ a) shows a view of the ecliptic plane, whereas Fig- 
ure[8]Jb) shows a viewpoint slightly tilted from solar west. 
Figure Etc) shows the viewpoint looking from the radial 
going through Earth toward the Sun. The orientations 
are visualized by cylinders with sizes derived from model- 
ing as taken from Table [3] However, in a manner similar 
to the correction we imposed for the shock surface, we 
tried to reproduce the sizes and positions of the MFR at 
one "snapshot" reference point in time, for which we de- 
liberately choose the S2 shock arrival time at the Wind 
spacecraft {ts2w)- To this end, we moved the flux rope 
centers according to Equation (1), but now using the av- 
erage speeds of the flux ropes (see Table |3J Column 5) to 
correct for the different arrival times at different heliocen- 
tric distances. As an example, the flux rope center arrives 
at VEX on August 3 01:30 UT, which is 15h35min ahead 
of the shock arrival time at Wind. During this time, the 
MS center moves 0.206 AU further away from the Sun 
(assuming a uniform speed of 550 km s" 1 for MS), and 
at the time ts2\v it is positioned at tm3v =0.933 AU, 
and this is the distance at which we place this MFR in 
Figure El For STEREO-B, the position of MS in Fig- 
ure |U (tmsv — 0.9381 AU) is less than re, because 
the center of the flux rope arrives later (August 4 02:26 
UT) than the reference time tsiw- In summary, at our 
reference time, the flux rope MS had an almost identi- 
cal heliocentric distance along the radials going from the 
Sun through the positions of STEREO-B and VEX. 

We also expanded (or contracted) the MFR size from 
the time of the observation of the MFR center at the 
spacecraft to the reference time, for all spacecraft ex- 
cept Wind. Similar to the estimate for the maximum 
magnetic field strength of the ICME flows at 1 AU as 
function of heliospheric longitude (see Equation ([3]) and 
Table 1), we employ a similar correction for the ongoing 
expansion of the MFRs. The average expansion of MCs 
in the inner heliosph ere, which we again use as a rough 
proxy, was quoted bv lLeitner et all (|2007l) as 

D{r) = 0.23 r 1 - 14 , (8) 

with D and r in (AU). A similar approach as in Section 
2.3 leads to 

D' = D h (r'/r h ) 1 - u . (9) 

This formula tells us the size to which an average MC 
would expand at any heliocentric distance r' if the diam- 
eter Dh of the MC is observed at a heliocentric distance 
r/j. As an example, for the application to VEX we use for 
r' the estimated position of the flux rope center (rM3v) 
from above and an estimated flux rope diameter at ry 
of Dy = 0.14 AU (again assuming the MFR has a uni- 
form speed of 550 km s _1 ). Equation ((9]) then leads to 
a diameter D' = 0.186 AU, which we use for the size 
of the pink cylinder at VEX in Figure [8] The result for 
STEREO-B is a contraction from 0.120 AU to 0.106 AU, 
which again follows from the fact that the center of MS 
hits STEREO-B later than t S2 w- 

4.4. Implications for Global Flux Rope Structure 
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TABLE 3 

Modeled Flux Rope Parameters at Each Spacecraft. 

FR Type Start Time At V(km s" 1 ) B z (nT) H m (°) m (°) D (AU) p (AU) * t (lO^ 1 Mx) $ p (lO^ 1 Mx) I z (M/ 

0) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14) 

STB M3 MV NES Aug 3 21:50 9h:12m 544 12.8 R -54 43 0.120 

VEX M3 MV -ES Aug 2 19:51 llh:21m ... 17.7 R -24 75 0.14 a 

Wind M2 GS SWN Aug 4 02:50 5h:08m 569 17.7 R 23 251 0.069 -0.004 0.071 0.383 228 

Wind M3 GS -WS Aug 4 09:38 13h:40m 544 16.6 L -39 301 0.157 0.032 0.127 0.570 47 

MES Ml MV NWS Aug 2 14:12 2h:29m 24.2 L 13 328 0.035 b .__ 



Note. — Column 1: spa cecraft, flux rope (FR), type of method: GS=Grad- Shafranov ( after IHu fc SonnenTp||2002t IMostl et al.M2009al) : MV 
— minimum variance (after Sonncrup & Cahill 1967: Bothmcr & Schwcnn 1998). Column 2: type of field rotation determined by eye (after 
IBothm cr & Schwcnn 1998; Mullig ari^t^dTrr998| : ) . Column 3: FR front boundary. Column 4: duration of the FR. Column 5: Norm of dcHoffmann- 
Teller velocity. Column 6: for GS: maximum value of the axial field component, MV: maximum magnetic field along the spacecraft trajectory. 
Column 7: FR handedness. Column 8: Inclination of the FR axis to the ecliptic plane (8 — [-90 ... +90]°). Column 9: FR axis longitude (cf>=[0 
... 360]°; sunward— 0°; solar cast— 90°). Column 10: GS: radial diameter D of the FR perpendicular to axis, MV: diameter along the spacecraft 
trajectory. Column 11: impact parameter p, the closest distance to the FR axis (only GS). Column 12: the axial flux of the FR (only GS). Column 
13: the poloidal flux for a FR length of 1 AU (only GS). Column 14: the axial current I z of the FR (only GS). 

a Assuming a radial speed of 550 km s" 1 times the durati on of the FR. For th e visualization of this FR size in Figure [8] an expansion to 
0.223 AU was estimated according to a quasi-linear law by Lcitncr et al. (2007) for the expansion of magnetic cloud diameters. The radial 
position of the FR in Figure [8] corresponds to the time of the shock arrival at Wind (see the text). 

b Assuming a radial speed of 600 km s — x . For the size of this FR in Figure[H] the Lcitncr et al. (2007) expansion law yields 0.075 AU at 
the time of the Wind shock arrival (see the text). 



In summary, Figure [8] should be seen as the final re- 
sult ol our flux rope analysis. We note several interesting 
hints on the global flux rope configuration, following from 
this figure. We indicated the MFRs which seem to be- 
long to the same large structure with similar colors: pink 
for MS at STEREO-B and VEX, blue for M2 at Wind, 
and red for Wind MS and MESSENGER Ml. These 
associations follow quite naturally by looking at the chi- 
rality of the field and the axial field directions. In par- 
ticular, MS for STEREO-B and VEX cannot extend to 
Earth because of the different axial field direction, which 
is eastward for the eastern (pink) structure and westward 
for the western (red) structure. Also, the field chirality 
is the opposite, so we think that these associations are 
solid. Further, it could be that Ml at MESSENGER 
belongs to the same structure as MS at Earth, because 
the axial field directions and chirality (left handed) are 
similar. However, these observations are very separated, 
in both space (0.6 AU radial distance, 49° heliospheric 
longitude) and time (about two days). While the sep- 
aration in longitude is less than an often quoted limit 
of 60° for the flux rope ext e nt in heliospheric l ongitude 
(IBothmer fc Schwennlfl998l) . I Webb et at] (120121) did not 
find connections between the two structures from SMEI 
density reconstructions. All this makes it difficult to tell 
if these observations sample the same flux rope or not. 

By looking at the longitude of the axes, best viewed in 
Figure [HKa), it seems that imagining the flux rope as a 
bent tube on a large scale works relatively well for these 
events. However, Figure [SJc) shows that the latitude of 
the axes for the pink and red tubes are quite inconsistent 
with each other - cautioning that the red tube might be 
a different structure. This might be a hint that an incli- 
nation of a flux rope should rather be viewed as a local 
parameter. It is, however, difficult to distinguish true ax- 
ial distortions or a "warped" shape on a large scale from 
a possible failure of the methods we used for finding the 
axial orientation in the first place. One should be espe- 
cially cautious of over-interpreting the orientation results 
if the field rotations inside the flux ropes are not well de- 
fined, as it is true here for the events at STEREO-B and 
VEX (see also Al-Had dad et al., 2012, S olar Physics, 
m revision ). However, IWebb et all (|2012D show fits of 



non-linear force-free models (jMulligan fc Russelll I2001D 
to these MFRs and they are quite consistent with the 
results presented here. Because different models give us 
quite similar answers here, we think that the picture pre- 
sented in Figure [8] is a relatively robust one. 

4.5. Associating ICME Signatures to CMEs 

In Figure [5] we try to relate the observed multi-point 
flux rope structures to the measured CMEs in corona- 
graphs. Similar to Figure Ha), orange arrows indicate 
initial directions for the four CMEs on 2010 August 1 as 
given by an elliptical forward modeling tec hnique applied 
to bo th STEREO /COR2 coronagraphs (jOdstrcil et al.l 
I2012T) . We think these results represent well the con- 
clusions of other papers to determine the direct ions of 
these CMEs. For example, iTemmer et al.l (|2012l ) found 
a range for the direc tion for CME2 from E6 to E41, and 
lOdstrcil et all (|2012l ) find E25 for CME2. Harrison et all 
(|2012l) . discussing the propagation phase in STEREO /HI 
for CME2 (in their paper called CME L), find a direction 
between E0 and E40, which is also consistent with E25. 
By using green ellipses, we group flux rope observations 
which likely can be traced back to a similar CME. 

First, it is quite straightforward to associate the MS 
structure at STEREO-B and VEX (pink tubes in Fig- 
ures M and H]) with the fast CME2, mainly because CME2 
was the one with the most eastward initial direction. 
It also had a significant northward direc tion of w N20 
(|Temmer et alj|2012t lOdstrcil et al.ll2012D . which should 
have contributed to a glancing encounter for these space- 
craft, resulting in the rather unpronounced magnetic field 
rotati on at STEREO-B and VEX. Note that IWebb et~aTl 
(|2012l) do not find a connection between these two struc- 
tures from SMEI density reconstructions, but it follows 
from the magnetic field modeling that they should belong 
to the same large-scale structure. 

It is intriguing that MS was seen at STEREO-B and 
VEX located 46° and 28° to the east of the most likely 
initial direction for the CM E2, respectively. Given the 
60° limit for widths of MCs (jBothmer fc Schwennlfl998l) 
noted earlier, one would expect VEX to encounter the 
edge of CME2 and that it misses STEREO-B completely. 
This flux rope is not observed at Earth, though it is 
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Fig. 6. — Plasma and magnetic field observations by the Wind spacecraft. From top to bottom: magnetic field strength (black) and GSE 
components Bx (red) and By (green), magnetic field B z (blue), proton bulk speed, proton density, proton temperature, proton /3, and 
the Dst (green line) and Kp (blue bars) indices for geomagnetic activity. The solid vertical line denotes the shock S2, and shaded intervals 
mark magnetic flux ropes. Two proton density enhancements are labeled as Dl/2. 
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Fig. 7. — (a, b) Transverse pressure relations and fitting func- 
tions (black solid and dashed lines), (c, d) Magnetic cloud cross- 
sections: the black contours are magnetic field lines in the x'-y' 
plane, which is perpendicular to the MFR axis, which is point- 
ing out of the paper plane. The spacecraft traverses from left 
to right, and the Su n is to the right (for exact definitions see 
Hu & Sonncrup 2002). The axial component of the field, B' z , is 
color coded. White arrows are magnetic field observations by Wind 
in the x'-y' plane, resampled to 15 points and used as initial val- 
ues for the numerical integration. Green arrows are small residual 
flow velocities in the deHoffmann-Teller frame. The white contour 
corresponds to the value of the vector potential A as indicated in 
panels (a, b) by a vertical line at At, and are defined as the bound- 
ary for the flux and current calculations (see Table |3j . 



25° west of the CME2 initial direction, which is closer 
than VEX. Explanations for this could be that CME2 
propagated more easily eastward into the depleted den- 
sity region by the 2010 July 30 CME, or that the in- 
teraction with CME1 around 30-50 R Q (|Liu et al.ll20T2l: 
iTemmer et aI|[20T2T) may have led to a deflection to the 
east for CME2 and to the west for CME1 (see also Lugaz 
et al. 2012, ApJ, in press, reaching similar conclusions 
for a different event). 

Indeed, at Earth, the M2 flux rope (blue), which does 
not have a counterpart elsewhere, might correspond to 
CME1, which had an initial direction also slightly to 
the east, but closer to Earth (see again Figure Fur- 
thermore, if the strong density peak D2 at the front 
boundary of MS is indeed filament material, then MS 
(the red tube) would correspon d to CME3, which was 
almost co-temporal w ith CME2 (jSchriiver fc Titlej|20ill: 
lHarrison et al.l |2012[ ) and associated with a filament 
eruption around disk ce n ter in the northern hemisphere 
(Schr iiver fc Title! [201lt ILi et al.l [MTU: lHarrison et all 
I2012D . As discussed in the previous section, while it 
makes sense to associate the only clean flux rope observed 
at MESSENGER with MS at Wind, the observations are 
too far separated in space and time to make a definitive 
connection. The weak CME4 did not seem to leave any 
strong trace in the in situ data, or, alternatively, could 



be responsible for the MESSENGER Ml flux rope. 

In summary, for three out of four CMEs we could find 
reasonable connections to MFRs observed in situ. For at 
least one CME, a deflection during heliospheric propa- 
gation seems to be needed to explain the connection be- 
tween the observations in coronagraphs and in situ. We 
have avoided associating the shocks with specific CMEs, 
because, for an extended shock surface like S2 spanning 
> 120° in heliolongitude, different CMEs may drive it at 
different frontal locations. 

5. A TWO-STEP GEOMAGNETIC STORM 

In this last section of the paper, we will take a closer 
look on this compound stream as observed near Earth, to 
find out which parts of it were responsible for producing 
geomagnetic activity. 

Figure [TU] plots from top to bottom (a) the clock angle 
of the interplanetary magnetic field (IMF), i.e., the polar 
angle in the GSM Y Z plane, (b) the Dst index at 1 bi- 
temporal resolution, and (c) the Kp index at 3 hr res- 
olution, obtained from the World Data Center for Geo- 
magnetism in Kyoto. The polar angle is defined as 0° for 
a purely northward pointing IMF, 90° for an eastward 
oriented, and 180° for a purely southward pointing field. 
The 1-hr Dst index is calculated from four low-latitude 
magnetogram stations and it monitors the strength of 
the symmetric equatorial ring current. In Figure UOf b) 
we have also plotted the pressure-corrected Dst index 
(Dst*) where the effect of the magnetopa use Chapman- 
Ferrar o currents is removed according to iBurton et al.l 
(|1975l) . This correction assumes a Chapman- Ferraro 
scaling and is not appropriate when the response of the 
magnetosphere to external forcing - i.e., to the interplan- 
etary motional electric fiel d (IEF) - becomes non-linear 
and approaches sa turation ([O'Brien fc McPherrorj|2000t 
iSiscoe et al.ll2002[ ). Saturation ty pically sets in at val- 
ues of the IEF ~ 6-8 mV m" 1 (|Hairston et al.l 120031: 
IMiihlbachler et al.l 120051) . In our case, IEF is less than 
6 mV m _1 most of the time (not shown) and thus a linear 
response may be assumed so that this correction should 
be valid. The Kp index is calculated from 13 subauro- 
ral stations and is thus more sensitive to high-latitude 
activity than Dst. 

Figure flOtb) shows that there are two main decreases 
in Dst*, with one significant recovery in between. The 
first depression was the strongest. The arrival of the 
shock and the following brief region of northward IMF 
and increased dynamic pressure were seen as a positive 
excursion in Dst*, which is known as a sudden impulse 
(SI). After August 3 19 UT, the IMF turned southward 
and Dst* started a rapid decrease. IMF was primar- 
ily negative throughout the sheath region and Dst* re- 
mained below -100 nT until August 4 03 UT after which 
it started to recover. The leading edge of M2 arrived 
at the magnetopause around August 4 04:30 UT, and 
thus the first Dst* depression was evidently caused by 
the turbulent sheath region, which, however, lead to a 
relatively steady southward pointing field, which could 
be a consequence o f the IMF draping around the ejecta 
in the sheath (e.g.. iGosling &; McComaj|1987t ILTu et al.l 
2008^). Th e Dst* minimum during this first depression 
was just below -100 nT at August 4 00 U T, thus it is clas- 
sified as an intense magnetic storm after iGonzalez et ID 
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Fig. 8. — Global configuration of magnetic flux ropes. Cylinders indicate sizes and orientations of large-scale magnetic flux ropes which 
could be successfully modeled. The red solid line is the approximate shape of the shock front S2. The positions and sizes of the MFRs are 
corrected to approximate the global configuration of the ICMEs at a snapshot in time, which we choose as the shock arrival time at the 
Wind spacecraft (2010 August 3 17:05). (a) View of the ecliptic plane from north, (b) Side view tilted to solar west, (c) Looking along 
the radial from the Sun to Earth. (An animation of this figure is available in the online edition of the Astrophysical Journal.) 
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(1994) The Kp maximum of 7 was reached on August 
3 21-24 UT. During the second peak the Kp maximum 
was 6+ from August 4 18-21 UT. 

During the passage of M2, the GSM B z component 
was mainly northward and Dst* recovered. A few hours 
after the arrival of MS at the magnetopause IMF turned 
southward and Dst* started its second decrease. Now, 
the IMF maintained steady southward orientation for 
more than 12 hr, keeping Dst* low. However, the mag- 
nitude of the southward component was only about 5 nT 
and no severe activity followed. Kp followed the same 
pattern as Dst*. 

As described above, the analyzed compound stream 
caused Dst* to develop in two distinct step s. Such two- 
step magnetic storms (K a mide et al.l Il998l ) are related 
to the double-structured southward IMF, caused typi- 
cally by a shock-driving south-north (SN) cloud, with 
the shock enhancing previously negative B z upstream 
fields. In such a case, the sheath fields initiate the storm 
that is further intensified by the f ollowing southwar d 
fields in the leading part of the cloud ([Vieira et al.| [2001). 
Shock-driving north-south (NS)-clouds may also lead to 
a two-step storm, but with a longer separation between 
two consecutive Dst* excursions and stronger recovery 
in between. If an NS-cloud has sufficiently long dura- 
tion, the ring current may decay during the first part of 
the cloud close to the quiet-time levels, and two sepa- 
rate m agnetic storms follow (according to Kamidc et al. 
(1998) this happens when the ratio of the magnitudes of 
the Dst* recovery and the first Dst* depression is larger 
than 0.9). However, the two-step storm analyzed in this 
paper was caused by the complex combination of a sheath 
region and two interacting MFRs. 

As the majority of flux ropes are expected to have NS 
field rotation during the rising phase of the current solar 
cycle 24 (e.g. 



19981 iLi et at 



Both 



mcr fc Schwennl ll998: Mul ligan et al.l 



201 If ) . it is interesting to ask how strong 



the ring current response would have been if M2 had 
NS-polarity instead of SN. To estimate thi s we used the 



modified Burton's Dst* prediction schem e (Burto n et al 



19751 iFenrich fc Luhmannl Il99l lO'Brien fc McPherron 
2000T ) for the interval of the compound stream with the 



sign of B z reversed during M2 (data not shown) . In this 
reversed -B z case, the compression by MS coincided with 
the southward fields in the end part of M2 and the fol- 
lowing storm had three distinct Dst* excursions, with 
the strongest depression occurring during the passage of 
M2. However, the Dst* minimum of the storm was only 
slightly stronger than when M2 had original polarity (-94 
and -70 nT for the reversed and real -B z M2 cases, re- 
spectively; note that the Dst* model underestimates the 
strength of the storms) and the magnitude of the third 
depression remained about the same (-64 nT and -58 nT, 
respectively). Although compound streams often cause 
intense Dst* storms (see Introduction) it is not clear how 
much pre-existing ring cu rrent affects the storm strength 
(e.g. jKozvra et al.ll2002t) . 

The studied event is also a good reminder of the im- 
portance of_th£_sheath fields i n generating magneti c 
storms (|Tsurutani et al.l 119881 : iHuttunen etall I2002T ). 
The strongest Dst* depression was caused by the sheath 
region, and not the magnetic flux ropes. The magnetic 
field magnitude in M2 was high as a result of interac- 
tion with MS, but the B z was mainly northward causing 



the recovery of the Dst* . Concerning MS, we already 
discussed that it is likely that the Earth made only a 
glancing encounter through this structure. As a conse- 
quence, the magnetic field magnitude was overall low and 
the following geomagnetic activity was not strong. 

6. CONCLUSION 

The series of events from 2010 July 30 to August 1 gives 
us a rare opportunity to "see beyond" single-spacecraft, 
in situ observations of ICMEs. We are actually unaware 
of any other series of events which was covered this well 
by four in situ spacecraft at almost equal longitudinal 
separations spanning 120° in longitude, which is a third 
of the heliosphere, though of course constrained to the 
ecliptic plane. As could be expected for a series of four 
coronal mass ejections heading clearly towards four out 
of five spacecraft from which we had the chance to obtain 
the in situ data, we found a plethora of physical phenom- 
ena such as shocks, flux ropes and interactions between 
various ejecta at these four locations, and attempted to 
link them into a global picture. 

We summarize our main findings as follows. 

1. A shock propagating partly into a previous ICME 
can lead to alterations of the s hock structure; a 
similar conclusion was reached bv lLiu et al.l ([2012). 
Additionally to this study, we have derived the 
approximate shape of the shock, showing it to 
be clearly non-spherical with respect to the Sun, 
though convex as expected. This non-spherical 
shape and a weakening of the shock at the eastern 
flank is caused (1) by parts of the shock propagat- 
ing into a medium carved out by a previous ICME, 
and (2) a higher initial CME speed of the eruptions 
on 2010 August 1 on the eastern solar hemisphere. 
From the shock normal orientations longitudes, the 
global shape can be reasonably extrapolated from 
single-spacecraft observations. This means that 
the longitude of the shock normal is pointing east 
(west) of the radial direction to the Sun for space- 
craft positioned on the eastern (western) flank of 
the shock, in three out of four shock normals (this 
was not true for the shock at MESSENGER). 

2. Concerning modeling of the magnetic flux ropes 
(MFRs), we found a relatively good consistency 
of the results fro m applying severa l methods (in- 
cluding those of iWebb et al.l 120121) to the MFRs 
at all in situ locations with the classical picture of 
ICMEs as large-scale, bent magnetic tubes (Fig. 
Et^a)). In particular, we were able to associate MS 
at STEREO-B and VEX as being parts of the same 
structure by looking at the chirality of the field and 
the axial field directions, at a longitudinal separa- 
tion of « 18°; in this way, it is also possible to 
associate MS at Earth with Ml at MESSENGER, 
though these observations are too separated from 
each other to make a definitive connection. While 
the longitudinal picture is consistent with the clas- 
sical paradigm, looking at the inconsistent inclina- 
tion of the flux ropes to the ecliptic plane (Figure 
[8jc)) might point to alterations of the flux rope pic- 
ture as "warped" rather than straight tubes. One 
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Associating CMEs to ICMEs 
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Fig. 9. — View of the ecliptic plane visualizing flux ropes as cy lindrical tubes and the shock surface S2 as a red solid line. The initial CME 
directions (orange arrows) are derived from STEREO /COR2 by Odstrcil ct al. (2012), with a shorter arrow indicating a later eruption. The 
STEREO/COH2 maximum field of view (FOV) of 15 Rq is shown by a yellow sphere. Green ellipses suggest which flux rope observations 
likely originate from the same CME. 
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may imagine holding a deformable "pool noodle"Q 
in one's hands - if deformed from its original cir- 
cular shape, the inclination to a given plane varies 
along the axis. We propose that we found some 
evidence for this axial "warping" in these events 
and thus, an inclination of a flux rope to the eclip- 
tic plane may be viewed ra ther as a local quant ity 
than a global one (see also iFarrugia et aT1l2011[) . 

3. A shock propagating into a previous eruption can 
give rise to a rotating magnetic field which can be 
perfectly modeled with various methods as a small 
cylindrical flux rope, but could be mistaken for a 
distinct flux rope originating from a CME {M2 at 
STEREO-B and VEX). 

4. This series of events do es not qualify as_a com- 
plex ejecta as defined by IBurlaga et al.l ([2002D be- 
cause individual flux rope structures can still be 
clearly identified and modeled at all spacecraft. 
However, some of the parameters have been sig- 
nificantly altered in the interactions that did oc- 
cur, i n accordance with IFarrugia fc Berdichevskv 
(2004). At Wind, in particular, a shock propa- 
gated through the preceding MFR, and/or ongoing 
energy and momentum transfer between the MFRs 
lead to heating and compression of the plasma (see 
also lLiu et al.ll2012D . 

In conclusion, we think that this series of events, as 
complex as they are at first glance, may qualify very 
well as a benchmark test case for space weather pre- 
diction models. In numerical simulations for forecast- 
ing CMEs, the outputs of the model runs should be 
consistent with the arrival times of the shocks and flux 
ropes, the number of MFRs at each spacecraft and the 
general profiles of plasma and magnetic field. First at- 
tempts in this direction using cone-shape d ejections were 
extre mely encouraging (jWu et all I201H lOdstrcil et all 
l2012f) , but need to include the flux rope structures in fu- 
ture simulations. This would make it possible to attempt 
predictions of the magnetic field components inside the 
ICME, and concerns in particular the ever elusive B z 
component (in GSM coordinates) which can be directly 
linked to the strengt h of the geomagnetic storm (e.g., 
iTsurutani et al.lll988f) . Situations like the one described 
in this paper may be expected to be more common during 
solar maximum. Sympathetic eruptions naturally lead to 
compound streams composed of several CMEs, altering 
ICME properties in distinct ways, sometimes pr oducing 
more geo-effective flows fe.g. JZhang et al.ll2007l ). 

New possibilities for serendipitous observations such as 
these are expected for 2015, when the STEREO space- 
craft will "cross" on the other side of the Sun as seen 
from Earth. With two additional spacecraft operating in 
the inner heliosphere at the moment, MESSENGER at 
Mercury and VEX at Venus, and the Sun then in its de- 
clining phase of the cycle, it might be possible to "catch" 
more ICME events on a global scale such as the one pre- 
sented in this paper. However, these future STEREO 
directed events might lack the three-dimensional imag- 
ing capabilites which can be used for the events in early 

For this term we would like to give credit to Volker Bothmer 
who first introduced it at a STEREO meeting in Dublin in 2010. 



2010 August, which is another reason why we think these 
events are a truly outstanding test case for the space 
weather community. 
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Fig. 10. — (a) Clock angle in the GSM Y Z plane of the interplanetary magnetic structure (0 ecliptic north, 180 ecliptic south), the (b) 
real time Dst index at 1 hr temporal resolution (black) and the corrected Dst* index (red), and the Kp index at 3 hr resolution. Vertical 
guidelines label shocks and magnetic structures. 
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